Direct measurement of finite-time disentanglement induced by a reservoir 



O 
O 

(N 

C 

a 

in 

(N 
> 

o 

(N 
0\ 
O 

in 
o 

a,: 
^— > 



X 



M. Franga Santos, 1 P. Milman, 2 L. Davidovich, 1 and N. Zagury 1 

1 Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Caixa Postal 68.528, 21.941-972, Rio de Janeiro, Brazil 
2 Laboratoire de Photophysique Moleculaire du CNRS, Univ. Paris Sud, 
Bdtiment 210-Campus d'Orsay, 91405 Orsay Cedex, France 
(Dated: February 1, 2008) 

We propose a method for directly probing the dynamics of disentanglement of an initial two-qubit 
entangled state, under the action of a reservoir. We show that it is possible to detect disentangle- 
ment, for experimentally realizable examples of decaying systems, through the measurement of a 
single observable, which is invariant throughout the decay. The systems under consideration may 
lead to either finite-time or asymptotic disentanglement. A general prescription for measuring this 
observable, which yields an operational meaning to entanglement measures, is proposed, and exem- 
plified for cavity quantum electrodynamics and trapped ions. 
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Entanglement is the most characteristic trait of quan- 
tum mechanics p| . As such, it has led not only to intense 
debate since the beginnings of quantum mechanics 0, 
but also to a variety ofpossible applications, ranging 
from communications 0, to computation [{| . Decoher- 
ence and entanglement are closely connected phenomena: 
not only decoherence follows from the entanglement of 
the system of interest with the rest of the Universe |(|, 
but also it is responsible for the fragility of entangle- 
ment in systems interacting with reservoirs. Because of 
this, understanding the basic mechanisms of decoherence 
and disentanglement has both fundamental and practi- 
cal implications. However, even for simple systems of two 
qubits, the available entanglement measures, like concur- 
rence as introduced by WoottersQ, involve mathemati- 
cal operations that do not seem to have a direct physical 
interpretation. It is desirable, therefore, to find exam- 
ples of systems for which an entanglement measure could 
be associated to an observable that could be easily mea- 
sured. This would not only make it easier to follow the 
dynamics of disentanglement for systems interacting with 
reservoirs, but it could be helpful for the physical inter- 
pretation of this process. Here we show that, for an ex- 
perimentally realizable example of a decaying system, it 
is possible to directly measure concurrence, through the 
detection of an observable that is invariant throughout 
the evolution of the system. 

A state of a bipartite system is disentangled or sepa- 
rable if it can be written in the form 



P = ^2 P i pi 



(1) 



For a pair of qubits, described by the density operator 
p, entanglement may be quantified by the concurrence^ 

C(p) =max{0, v 7 ^- V^- - \/M , (2) 

where the A,'s are the eigenvalues, in decreasing order, 
of the Hermitian matrix p(a y <g> o y )p*{<j y ®a y ). It ranges 



from C = for a separable state to C = 1 for a maximally 
entangled state. 

The connection between disentanglement of bipartite 
systems and decoherence has been recently analyzed, 
within the framework of models that may lead to finite- 
time disentanglement, while the individual subsystems 
decay asymptotically in timejil 1^. fiol HH . These ex- 
amples are interesting insofar as they allow a clear separa- 
tion between the disentangling and the decay dynamics. 
In [lfl j. a realistic system was considered, consisting of 
a pair of two-level atoms, each interacting with its own 
reservoir. Depending on the initial atomic state, one may 
have cither finite-time or asymptotic disentanglement. 

Simple methods to detect this transition from entan- 
glement to disentanglement would be clearly desirable. 
However, the non-physical nature of the operations in- 
volved in the definition of concurrence imply that such 
a detection is usually a challenging problem, requiring 
the full tomographic reconstruction of the state. Entan- 
gled states may also be identified through entanglement 
witnesses, which are non-positive operators that are pos- 
itive in the subspace of separable states 01 . Entangle- 
ment witnesses allow one to identify some, but not all, 
entangled states. 

For two trapped ions, entanglement between corre- 
sponding two-level internal states has been established 
either by means of a witness involving a sin gle matrix 
element of the two-qubit density operator [lil ll.'jj , or by 
measuring the full concurrence via tomographic recon- 
struction of the state ^(| . 

In this paper, we show that it is actually possible to 
detect, through the measurement of a single observable, 
finite-time disentanglement of an initial pure state of 
two qubits, which evolves under the action of a reser- 
voir. This observable is a "perfect witness" for the 
class of states here considered: any entangled state in 
this class would lead to a negative value for this opera- 
tor, which is precisely equal to minus the concurrence 
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of the state. Remarkably, the same observable yields 
the concurrence throughout the evolution of the sys- 
tem. It is possible therefore, by measuring this operator 
as the system evolves, to pinpoint the precise moment 
when disentanglement occurs. Also, a simple physical 
interpretation can be given to concurrence in this case. 
Our proposal is within present experimental capabili- 
ties in systems like trapped ion s 1 14 
electrodynamics (cavity QED)|l71 
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electrodynamics [l^f, and nuclear magnetic resonance|20|. 

We consider here initial states of the form |^(0)) = 
H|00) + |/3| exp(i6>)|ll), where and 1 correspond to 
the ground and excited state of each qubit, respectively. 
The two qubits may stand for spin-down and up internal 
states of two trapped ions, or to one- and zero-photon 
states of two modes of the electromagnetic field in the 
same or in different cavities. States like this have been 
realized in experiments with trapped ions [l^lllfil Ell l2^ | . 
The phase 9 is attributed by the preparation process. We 
assume that the two qubits evolve under the influence of 
low-temperature independent and identical reservoirs, so 
that higher-energy states are not populated. 

Under these conditions, we show that the concurrence 
of the density operator evolving from the above initial 
state is equal to C(t) = max{0, — W(t)}, where W(t) = 

Tr 



p{i)Wg , and We = 1- 2|$(0))($(0)|, with 
= (|00)+e 4e |ll))/\/2. 



(3) 



The observable Wg is a "perfect witness," the same 
throughout the evolution of the system: it is positive for 
separable states, and negative for entangled states. 

We note that W{t) = l-2P(i), where P(t) is the prob- 
ability of finding the system in the state at time 
t. It is easy to show that other choices of the relative 
phase in (J2J, different from 9, would yield smaller val- 
ues of the corresponding probability. This yields a sim- 
ple physical interpretation of concurrence in this case: it 
is twice the maximal excess probability, with respect to 
50%, of finding the state of the system in a maximally 
entangled state corresponding to the subspace spanned 
by {|00), |11)}. One should note that evaluation of W(t) 
amounts to measuring this probability. We show now 
that this can be done in a simple way, by inverting the 
process which yields |$(0)) from the state |00). 

State may be obtained from the state 1 00) by ap- 
plying to it an unitary transformation U (9) composed of 
a one-qubit ir/2 rotation followed by a CNOT operation. 
Therefore, the probability P(t) can be written as 

P(t) - (<!>(9)\p(t)M0)) = (00\U~ 1 (9) P U(9)\00) , (4) 

that is, the probability of finding the system in the state 
can be obtained by applying to the system the inverse 
of U(9), and then detecting the probability of finding 
both qubits in the ground state. 



The general format, in the basis |11), |10), |01), |00), 
of the time-dependent density matrix that evolves from 
state |*(0)) is: 



(w{t) z{t)\ 

x{t) 

x(t) 

\z*(t) y(t)J 



(5) 



with w(t), x(t), and y(t) real. Indeed, the incoherent 
decay of the state |11) under the action of the reservoir 
does not lead to coherence between the states |10) and 
1 01) . Also, the symmetry of the initial state and the 
equality of the damping rates for both reservoirs implies 
that the population of theses two states is always the 
same. 

The concurrence corresponding to this state is easily 
determined as C(t) = max{0, 2\z(t)\ - 2x(t)}. There- 
fore, the state is entangled if and only if \z(t)\ > x(t). 

The probability of finding the system described by (JjJJ 
in the state J2J is given by 



P(t) = [1 - 2x(t)+2\z(t)\]/2, 



(6) 



and therefore C(t) = max{0, 2P(i) — 1}, as anticipated. 

For the initial state |*(0)>, C(0) = 2\a/3*\ = 2\z(0)\, 
which equals one for the maximally entangled state 
(| 00) + e i6 \U))/y/2. As time evolves, the state becomes 
disentangled when P(t) reaches the value 1/2. 

This result is rather insensitive to small deviations 
from the above state. Thus, if the populations of states 
1 10) and 1 01) are slightly different (due, for instance, to 
unequal decay rates), then if e(t) = \xi(t)—X2(t)\ -C \z{t)\ 
it is easy to show that the concurrence is still given by 
the above expression, up to terms of order e 2 . Also, our 
method can be extended to any state in the subspace 
{|00>,|11>}. 

We show now that, for the usual linear coupling with 
markofnan reservoirs, C(t) may reach zero at finite times. 
The master equation governing the time behavior of the 
system can be written in the Lindblad form |2^ |: 



/5c -c, 



(7) 



where 7j are decay rates and c, , cj are operators repre- 
senting the coupling to the reservoir. For atoms decaying 
spontaneously due to the coupling to a zero temperature 
reservoir, c = <7_, where <r_ is the lowering operator act- 
ing on the electronic levels of the atom. For modes of 
the electromagnetic field in one or two cavities, coupled 
to a zero- temperature reservoir, the decaying dynamics 
is analogous to the ionic system, provided that the pho- 
tonic state is initially in the subspace {|0) , |1)}. In this 
case, Cj correspond to annihilation operators for photons, 
for each of the two cavity modes. The solution of (0 for 
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FIG. 1: P as a function of r = ft for four different initial 
states, of the form a|00) + /3\11), with concurrence C(0) = 
2|a||/3|. Solid thick lines correspond to C(0) = v2/3, while 
dotted lines correspond to C(0) = VT5/16. Concurrence, as 
a function of time, is given by max {0, 2P(t) — 1}. Entangle- 
ment persists as long as P > 1/2 (this bound corresponds 
to the thin solid line). States with the same initial entan- 
glement exhibit finite-time disentanglement for \/3\ > \a\ and 
infinite-time disentanglement for \a\ > \0\ 

any initial state lying in the subspace {|00) , |H}} is of 
the form JSJ with: 

w(t) = w(0)e- 2 ^, 

x(t) = w(0)(l-e- 7 ')e- 7 *, (8) 
y(t)=y(0) + w(0)(l-e-^) 2 , 
z(t) = z(0)e^ 7 *. 

Consequently, an initial state of the form a 1 00) + j3 |11) 
becomes separable for t s = — ~ln(l — fej). It is clear 
that only for states with > \a\ is this condition ful- 
filled for finite times. The probability P(t) is shown in 
Fig. H for four different initial states. One should note 
that states with the same initial entanglement exhibit 
finite-time disentanglement for |/3| > \a\ and infinite- 
time disentanglement for \a\ > |/J|. For states of the 
form 7|01) + <$|10), disentanglement occurs asymptoti- 
cally in time. This justifies our interest in the special 
class of states here considered, since they allow a clear 
separation between the processes of disentanglement and 
decohcrencc, besides allowing a simple monitoring proce- 
dure of the disentanglement process. 

One should note that the method of measurement out- 
lined above is valid only if the operations on the qubits do 
not involve auxiliary systems, as is the case for instance 
in nuclear magnetic resonance. On the other hand, the 
realization of quantum gates in ion traps (where the vi- 
brational mode is used for the realization of a CNOT 
gate) , or in cavity QED (where the cavity mode mediates 
the interaction between two atoms, or an atom mediates 
the interaction between two cavity modes) require auxil- 



iary systems. In these cases, the operation U(9) does not 
depend only on the qubit operators, so that Eq. (@J does 
not hold. In spite of this, the above method can still be 
used to measure entanglement, the results differing from 
the concurrence by a scaling constant. We exemplify this 
with application to ion traps and cavity QED. 

For two trapped ions, the production of the initial state 
and the measurement of P(t) can be done by employing 
standard techniques involved in t he p roduction and de- 
tection of entangled states [la, |2j|, (22j . 

In these experiments, maximally entangled states are 
produced through a combination of red-sideband and 
blue-sideband transitions and detected using the shelving 
technique: the ion undergoes cycling transitions to an ex- 
cited unstable state if and only if it is in the ground state, 
the resulting fluorescence being measured with very high 
detection efficiency (of the order of 99%). The fluores- 
cence yield of each ion is thus proportional to the prob- 
ability of finding it in the ground state. 

One should note that the decay of the ions does not 
lead to population of the vibrational modes, which re- 
main in the ground state. We consider the initial state 
(|a||gig20) + e l6 \{3\\eie 2 0)) / where and gi stand re- 
spectively for the upper and lower internal states of the 
ions (i = 1, 2) and integer numbers stand for the states of 
the vibrational mode. Disentanglement of this state can 
be measured by applying to the first ion a red-sideband 
pulse, so that |gil) — > |ei0) and |ei0) — ► — | (71 1) . This 
operation transforms the state p{t) ® |0) (0|, with p{t) 
given by (0, into 

w(t) \gie 2 l) (gie 2 l\ + y(t) \gig 2 Q) (5152OI 
-z(t) \ 9l e 2 l) (gm0\ - z*(t) \ gi g 2 0) ( 9l e 2 l\ 
+x(t){\ gi e 2 0) ( 9l e 2 0\ + \g l92 l) (gig 2 l\) . 

Next a bluc-sideband pulse is applied on the sec- 
ond ion, so that \g 2 0) -> (| 52 0) - exp(i<S)|e 2 l))A/2, 
leal) (exp(i5) | fl2 0) + |e 2 l))/V2, and \g 2 l) 

cos(?rV2/4) \g 2 l) -sin(7r^/4) exp(-iS) |e 2 2). The prob- 
ability P gg ((5, t) that both ions are in the ground state is 

P gg (S, t) = 1/2 - \z(t)\ cos(# -5)- rjx{t) , (9) 

where r\ = sin 2 (7T-\/2/4). Choosing S so that cos(# — S) = 
—rj one gets, calling V gg {t) the value of P gg (S,t) at this 
point: 

2V gg {t)-l = r 1 [2P{t)-l], (10) 

with P{t) given by I©- This shows that max{0, 2V gg — l} 
is proportional to the concurrence. One should note that 
even if the initial phase 9 is not known, measurement 
of P gg (d,t) for three values of the phase S, leading to 
three linearly independent equations for x(t), \z(t)\, and 
9, would allow one to determine the concurrence. 
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FIG. 2: Experimental setup for probing finite-time decoher- 
ence in cavity QED. 



For cavity QED, the experimental setup involves two 
high-Q cavities (C a and Cb) on either side of a low-Q cav- 
ity (Caux), as shown in Fig. [21 A two- level atom crosses 
this system, interacting resonantly with the three cavi- 
ties. The field in C aux can be taken as classical, gener- 
ating a 7T rotation of the atomic state. The fields in C a 
and Cb are initially in the vacuum state. The interaction 
time between the atom and each cavity is adjusted, for 
open-cavity geometry [tR flij]. by Stark-shifting the rel- 
evant atomic levels, so as to tune them into resonance 
with the cavity mode for the proper amount of time. 

The first step is the preparation of the entangled two- 
cavity state. This is done by sending a two-level atom, 
initially in the excited state |e), through the three cav- 
ities. The interaction time with the first cavity (C n ) 
is adjusted so that state |<Fi) = a\e, 0, 0) — (3\g, 1,0) 
is produced (\i,n,m) denotes atomic state \i) and Fock 
states \n,m) for cavities C a and Cb respectively). That 
is, the atom has a probability \a\ 2 of remaining in the 
initial state, and a probability |/3| 2 of decaying to the 
state \g), leaving one photon in the resonant mode of 
C a . The atomic 7r rotation in C aux , takes into 
1*2) = ot\g, 0, 0) + /3|e, 1, 0), and another ir rotation, this 
time with the resonant mode of Cb, produces the desired 
entangled state /) = a\0, 0) +/3|1, 1) of the two modes. 
The atom leaves the setup always in the ground state. 

Due to cavity losses, the prepared state f) evolves 
into the mixed state described by Eq. ©. In order to 
measure the disentanglement of the state of the high-Q 
cavities, another two-level atom is sent again through the 
setup. This atom enters C a in the ground state and un- 
dergoes a Rabi 7r rotation if there is one photon in the res- 
onant mode. Next, C aux is used to generate yet another 
atomic 7r rotation and then the atom state undergoes a 
final rotation in Cb, so that \gl) — ► (\gl) + |e0))/v2, 
l e °) (- Iff 1 ) + |e0))/\/2. The rotation in C aux includes 
a control phase 5/2 (\g) -> e~ 4,5 / 2 |e) and |e) -> -e lS / 2 \g)). 
Finally, the atomic internal state is measured by ioniza- 
tion. The probability P e (6,t) of finding the atom in the 
excited state is given by the same expression obtained 
for the two-ion system, leading to the concurrence in the 
same way. 

Finite-time separability is related, in the example dis- 
cussed above, to properties of the reservoir and the initial 



state. The same phenomenon occurs for a diffusive reser- 
voir, acting on the two-qubit system. The corresponding 
contribution for the master equation can be written in 
the form J7J) with c.\ = rj_ and c 2 = rj + for each reser- 
voir. Due to the symmetry of this reservoir, all states be- 
longing to the subspace {|00) , |11)}, or yet {|01) , |10)}, 
display finite-time separability, independently of the rel- 
ative population of the states. This is consistent with the 
findings in which considered the action of classical 
noise on a two-qubit system. 

Since the work of John Bell [3| the subtle property 
of entanglement has been subjected to many quantita- 
tive tests and has led to intriguing consequences. Under- 
standing the physical meaning of entanglement measures 
remains however a major challenge. In this paper, we 
have analyzed an example of a decaying two-qubit sys- 
tem for which it is possible to attribute an operational 
meaning to an entanglement measure, valid throughout 
the decay process, and we have proposed an experimen- 
tal procedure which amounts to a direct measurement of 
entanglement. 

We acknowledge support of the Millennium Institute 
for Quantum Information, FAPERJ and CNPq. 
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